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Abstract 

The aim of this paper is to summarize some recently obtained relations between the 
Ablowitz-Ladik hierarchy (ALH) and other integrable equations. It has been shown that 
solutions of finite subsystems of the ALH can be used to derive a wide range of solutions 
for, e.g., the 2D Toda lattice, nonlinear Schrodinger, Davey-Stewartson, Kadomtsev- 
Petviashvili (KP) and some other equations. Similar approach has been used to construct 
new integrable models: 0(3,1) and multi-field sigma models. Such 'universality' of the 
ALH becomes more transparent in the framework of the Hirota's bilinear method. The 
ALH, which is usually considered as an infinite set of differential-difference equations, 
has been presented as a finite system of functional-difference equations, which can be 
viewed as a generalization of the famous bilinear identities for the KP tau- functions. 

1 Introduction 

One of the characteristic features of the theory of integrable systems is the fact that various models, some- 
times rather different apparently, turn out to be, in such or other way, closely interrelated. For example, the 
nonlinear Schrodinger equation (NLSE) is gauge equivalent to the Landau-Lifshitz equation, sine-Gordon 
to Thirring model, Ablowitz-Ladik model - to the classical Heisenberg chain, etc. A large number of the 
integrable equations can be obtained as reductions of, e.g., the KP or the self-dual Yang-Mills equations. 
Poles of the rational solutions of the KdV evolve according to the Calogero model and so on. 

Another kind of interrelations has been discovered by Levi, Benguria jl], ||, Flaschka ||], and also Shabat, 
Yamilov They considered sequences of Backlund transformations (BT) for some integrable nonlinear 
problems 

iV[V'] = 0, ... ^ V ^ V-i ^ ^ - 

and demonstrated that these sequences can be described by differential-difference equations (DDE) which 
are also integrable. For example, sequences of BT's for the NLSE provide solutions for the Toda chain. 

Similar ideas, in a somewhat transformed form, have been used in the works 0. It turns out that 
in some situations solutions of the DDE's can be used to generate families of solutions for some related 
partial differential equations (PDE). Namely such cases have been discussed in l^j, |[ 0]. There the approach 
of H, Q has been extended by taking into account several DDE's simultaneously, which enabled to deal 
with multidimensional PDE's, such as the 2d Toda lattice (2DTL), and the Davey-Stewartson (DS) equation 
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(2+1-dimensional systems) and others. The method of the works j|, ^, Q can be explained by discussing the 
following very simple examples. 

Toy example 1. 

Consider the system of differential-difference equations 

idxQn = Qn+l + q-n-l -|^^ 

By trivial algebra one can show that, first, this system is compatible and, second, that for every n the 
quantity q = Qn satisfies the Helmholtz equation. 

Aq + 4g = 0, A^dxx + dyy (1.2) 

This example is aimed to demonstrate that the system of DDE's can be converted to PDE. In the linear 
case this is not very interesting, but in the nonlinear one, as will be shown below, such transformations may 
be fruitful. 

Toy example 2. 



As the second example consider the system (1.1) enlarged with one equation more: 



idxqn = Qn+l + Qn-l 

dyqn ^ q-n+l - q-n-l (1-3) 
idtqn = qn+2 + gri-2 

Again, one can straightforwardly verify that this system is compatible and obtain that for every n the 
quantity q — qn satisfies the equation 



idtq+-Dq = 0, a^dxx~dyy (1.4) 



1 

-I 
2 

As in the first example we have converted the system of DDE's into a PDE. But it should be noted that we 
started with three 1+1, i.e. 2 dimensional equations, while the resulting one is 2+1, i.e. 3 dimensional. This 
is not crucial in the linear case, but is very important in the nonlinear one. The methods elaborated for mul- 
tidimensional {d > 3) nonlinear integrable systems (2DTL, DS, KP, etc.) such as, e.g., the multidimensional 
inverse scattering transform (MIST), are much more complicated than the traditional one (1ST) developed 
for the 2- and (l+l)-dimensional systems (KdV, NLSE, Ablowitz-Ladik model, etc.). It will be shown below 
that in the nonlinear case such 'decomposing' of a multidimensional PDE in a system of low-dimensional 
DDE's can be rather useful, at least from the practical viewpoint. 

Transformations similar to the ones described above, their nonlinear variant, are the core of the works 
1^, ^, . But before proceeding further I would like to discuss the following very important question, namely 
the question of compatibility. It is rather easy to write down compatible linear systems as the ones above, 
but in the nonlinear case the problem is a little bit more difficult. Nevertheless it is surely possible to 
construct systems, say, ones similar to (p"!]), (p~3|). But there is no need to do that, because such nonlinear 
and compatible systems are already known ~ they are the integrable hierarchies. Every integrable equation 
does not appear alone, it is always a member of some infinite set of related equations. And commutativity 
of corresponding flows, or involutivity of corresponding Hamiltonians, or, in simpler words, compatibility of 
corresponding equations, is one of the ingredients of integrability. Thus, in what follows I will 'nonlinearize' 
the toy examples using an integrable hierarchy, namely the Ablowitz-Ladik hierarchy. 

The plan of the present paper is as follows. After outlining some basic facts related to the ALH (section 
H) I will establish relations between this hierarchy and some other equations, such as 2DTL, DS, NLSE, KP 
(section^) as well as few vector models (section Q). To make this 'universality' of the ALH more transparent 
I will discuss it in the framework of the Hirota's bilinear method and present in section ^ the recently 
obtained § functional representation of the ALH which can be viewed as a generalization of the famous 
bilinear identities for the KP tau-functions. 
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2 Ablowitz-Ladik hierarchy. 



The ALH is an infinite set of ordinal differential-difference equations, that has been introduced by Ablowitz 
and Ladik in 1975 The most well known of these equations is the discrete nonlinear Schrodinger equation 
(DNLSE) 

iqn = qn+l - '2qn + Qn-l " f/n^n (^n+l + Qn-l) (2.1) 

and the discrete modified KdV equation (DMKdV) (see, e.g., 

Qn = Pn {qn + l - qn-l) (2.2) 

where 

p„ = 1 - (?„r„, r„ = -KQn, K = ±1 (2.3) 
All equations of the ALH can be presented as the compatibility condition for the linear system 

= Un^n (2.4) 

dt'fn = (2.5) 

where dt stands for d/dt, which leads to their zero-curvature representation: 

dtU,,^Vn+lUn-U,,Vn (2.6) 

In the standard 1ST approach developed in Q the matrix L/„ for the ALH is given by 

where A is the auxiliary (spectral) constant parameter. 

According to j9|, elements of the matrix Vn, can be chosen as Laurent polynomials in A in such a way that 



(2.6) holds automatically for all A's provided g„'s and r„'s satisfy some differential relations. It should be 
noted that one can obtain an infinite number of the matrices Vn (which are Laurent polynomials of different 
order) which leads to the infinite number of differential equations dqn/dt = F^, {I = 1,2,...). Using the 
widely accepted viewpoint one can consider g„'s and r„'s as depending on the infinite number of 'times', 
qn = qn{ti,t2, ■■■) and consider the Ith equation of the ALH as describing the flow with respect to the 
Ith variable, dqn/dti = F^. Traditionally it is implied that all 'times' ti are real, which is grounded from 
the standpoint of physical applications, and also is convenient in the framework of the inverse scattering 
technique. However, in some cases it is more convenient to use instead of real 'times' U some complex 
variables Zj,Zj, j = 1,2,... (as in, say, 2D Toda theory), which, as will be shown below, exhibit in a 
more explicit way some intrinsic properties of the ALH. A simple analysis yields that the family of possible 



solutions of (2.6) (and hence the equations of the hierarchy) can be divided in two subsystems. One of them 
consists of V^-matrices which are polynomials in A~^ (I will term the corresponding equations as a 'positive' 
part of hierarchy), and the other consists of matrices which are polynomials in A ('negative' subhierarchy) , 
while in the standard, 'real-time', approach all V-matrices contain terms proportional to A™ together with 
the terms proportional to A"™ (m > 0). Let us consider first the 'positive' case. An infinite number of 
polynomial in 1/A solutions {j — 1,2, ...) possesses the following structure: 

In n 



where the elements al^, (5^ satisfy the equations 



ai+i - < = -qnPi+i + Tnli (2-9) 

- Si = qnPi - rnlUi (2-10) 

djqn = qnSi+i + 7;^+! = qnai^^ + 7ri^^ (2-11) 

djrn = -TnSi ^f3i = -r,,aiX\ - PiX\ (2-12) 
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with dj = d/dzj. Choosing 
we can obtain consequently 



= 6° = = 0, d° = -z (2.13) 



a,\ =0 al = -irn-iqn 



(2.14) 



and, in principle, all other matrices V^. This leads to the infinite system of equations for q„, r„, some first 
of which are 

diq„ = -ipnqn+i (2-15) 

diTn = ir„_ip„ (2.16) 

d2q7i = -WnPn+iqn+2 + irn-lPnqnqn+l + iPnTnqn+l (2-17) 

d2rn = irn~2Pn^lPn ~ «f n-lPn?'ng«+l - «r^_iP„g„ (2-18) 

Analogously, looking for the IZ-matrices of the form 

Kr^ = A^Kr^- + (;^ ^) (2.19) 

one can obtain the 'negative' part of the ALH. Some first of its equations are 

diqn = -iq„-ip„ (2.20) 

9ir„ = ipnrn+i (2-21) 

d2qn = -iqn-2Pn-lPn + iqn-lPnqnrn+l + iql^iPnTn (2.22) 

d2rn = ipnPn+irn+2 - ig„-ip„r„r„+i - ipuqurl+i (2.23) 

where dj — d-j — d/dzj and Zj is the complex conjugated of Zj. 

I will not discuss here this hierarchy in detail, because now we have all the necessary to derive relations 

between the ALH and other integrable models and namely this is the main theme of the present paper. 



3 'Embedding' into the ALH. 



3.1 The Ablowitz-Ladik hierarchy and the 0(3,1) cr-model. 

The first, and the simplest, implementation of the 'embedding' into the ALH method was carried out in |5| 
and can be viewed as a nonlinearized version of the first toy example. 
Consider the following system of two equations of the ALH: 

ida:qn = Pn (^ri+l + qn-l) /o 
dyqn = Pn {qn+l ~ qn-l) 



The first of them is the DN LSE (2.1) transformed by means of the s ubstit ut ion q n qn exp{2ix), while the 



second is the DMKdV (2^). One can also view these equations as ( ^.15 ), ( ^.20[) rewritten in terms of the 



variables x = Re zi and y = Im zi . 

This system is compatible because it is a system of equations belonging to one hierarchy (this fact is not 
hard to verify directly in this case). Differentiating the first equation with respect to x and the second one 
with respect to y one can get the identity 

div — Vg„ + 2 (p„_i +p„+i)9„ = (3.2) 

Pn 
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from which, using again (3.1), one can derive the foUowing one: 



1 



div — Vg„ 



Thus, we have obtained that for all n's the quantities q — qn solve the PDE 



1 



4 - ^ |Vg„| 

Pn 



Aq - K- (Vg, Vq) + ipq = 
P 



where 



p= l + K\q\ 



(3.3) 



(3.4) 



(3.5) 



which is a nonlinear analog of the Hclmholtz equation (1.2). It turns out that this equation is not merely 
some 'nonlinearization' of the Helmholtz equation. It has a rather interesting physical origin. This is the 
field equation of some cr-model which has been proposed and discussed in |^ and which should be called, 
using the terminology by Pohlmeyer [Ol, the 0(3,1) cr-model. This model arises from the problem 



under the restriction 



Here A is the two-dimensional Laplacian, 



fi'^ is a space-like vector from the Minkowski space, 

AfJ- — ( Afi ^2 /3^ 



with the scalar product 



3„;,3 



>^ - (/-^V - 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



To satisfy the condition ( |3.7[ ) the Lagrange multiplier A is to be set to A = — (V0^, V(/)^). The vectors 
(/)^, dxcf)^, dy4>^ together with the time-like unit vector which is normal to the surface (j)^{x,y), form a 
local basis in the Minkowski space which satisfies the Gauss- Weingarten system. 



(3.10) 



where IJ and V are some 4 x 4-matrices (not written here). The integrability conditions for the system 
( 3.1C ), the so-called Gauss-Codazzi equations, can be presented as follows (see [|| for details): 



and 



COS (X 

Aa - 4 sin a cos a ( V/3, V/3) = 

sin a 



div (cot^ a V/3) = 



(3.11) 



(3.12) 



(note that the particular case of the above equations, /3 = 0, is the well-known elliptic sine-Gordon equation). 
These equations are the Euler-Lagrange equations for the action 



dxdy C 



with the Lagrangian 

C = (Va, Va) + cot^ a (V/3, V/3) - 4 cos^ a 
which can be rewritten in terms of the function 

q = cos a ■ exp(i/3) 



(3.13) 

(3.14) 
(3.15) 



as 
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This Lagrangian can also be obtained as the reduction of the Euchdean version of the principal chiral 
field model, as is outlined in the Appendix of 

The field equation corresponding to the Lagrangian ( 3.1(^ ), 



Ag+Y-^(Vg,V<7)+4(l-|g|^)g = (3.17) 



is nothing other than equat ion (3.4) for k = —1. 

The field equation ( 3.17 ) is integrable and it is possible to develop the inverse scattering scheme applicable 
to this equation (it has been done in the paper ||l^ ). However, it turns out that the derived relationship 
between this model and the ALH provides almost all results that are usually obtained in the framework of 
the ISM. It can be used to demonstrate that the 0(3, 1) cr-model possesses an infinite number of symmetries 
and conserved quantities and to obtain soliton and some other solutions for the field equation. 

3.2 The Ablowitz-Ladik hierarchy and the 2DTL. 

In this section I want to present some results that have no nontrivial linear analogues. Let us consider agai n 



the D NLSE-DMKdV system, now written in terms of the complex variables zi and zi, i.e. the system ( 2.15 ), 



(|2.20| ) 



= -iPnqn-l 



(here d stands for d/dzi and d for d/dzi) and turn our attention to the quantity p„. One can derive from 

(3.19) 



(3.18) the corresponding equations for p„ 



dpn = inpn (qnqn+l - qn-iqn) 

and obtain from this system, together with ( ^.1§| ), that the quantities Pn satisfy the following equation: 

ddXnpn = Pn-i - 2p„ + Pn+l (3.20) 

This equation is the famous 2DTL equation which can be rewritten in terms of the functions u„ defined by 
Pn = exp (u„ — Un+i) and the real variables x and y (zi = x + iy) as follows: 

^Alin = exp (li„_i - Un) - exp (m„ - Un+l) (3.21) 

Thus, we have shown that the 2D Toda lattice turns out to be hidden in the simplest equations of the ALH 
or, in other terms, the 2DTL can he 'embedded' into the ALH. 

The main, so to say, 'practical' value of the ALH-2DTL correspondence is in the fact that we can obtain 
a wide number of results for the 2DTL, which is 2+1 dimensional system, by means of the traditional 
version of the inverse scattering transform which has been developed for the 1+1 dimensional systems such 
as DNLSE, DMKdV, without invoking the MIST (that has been elaborated for the 2DTL by Lipovsky V.D. 
and Shirokov A.V. in Q). AU known solutions for the ALH (read for the DNLSE or DMKdV) can be 
converted to solutions for the 2DTL. 

Another interesting problem is the question of the conservation laws jT^ . One can now obtain an infinite 
number of the conservation laws in terms of the ALH and then 'convert' them to the 2DTL case (this idea 
has been described in Q). Some new (comparing with the paper p5[|) conservation laws have been obtained 
in Q using the following simple procedure. Take one of them 

divJ = (3.22) 

and apply operator d/dtj where tj is one of the hierarchy's times. This will yield an infinite number of 
conservation laws 

div J J = (3.23) 
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An interesting manifestation of the ALH-2DTL relation is the fact that the so-cahed conserved quantities 
'of discrete variable direction' which are of the form 



dl = 



(3.24) 



derived by Kajiwara and Satsuma 1 15 when considered in the framework of the approach of are nothing 
other than constants of motion of the equations of the ALH, i.e. Hamiltonians of the ALH-flows. 



3.3 The Ablowitz-Ladik hierarchy and the DS equation. 

Now let us discuss the bigger subsystems of the ALH, and start with the nonlinearization of the second toy 
example. Consider the system (2.15) - (2.22) 



diqn = ~ipnq,i+i 

diq„ = -iqn-iPn 

d2qn = -iPnPn+iqn+2 + ir„_ ip„g„(7„+i + ipnrnql+1 

d2qn = -iqn-2Pn-lPn + lPn'7n?'n+l + iql-lPnTn 



Differentiating ( p.25| ) with respect to zi one can straightforwardly obtain, using ( 3.27 ), that 

id2qn + dlqa = -2C„q„ 



where C„ — r„_ip„g„+i. On the other hand, it follows from ( 3.26| ) that 

diCn = iPn {rnqn+l ~ r^-iqu) = diPn 



which leads to 



(3.25) 
(3.26) 

(3.27) 
(3.28) 



Here I do not write the corresponding equations for r„'s, which can be restored from ( 3.25 ) - ( 3.28| ) using 
the involution r„ = — 



(3.29) 

(3.30) 

(3.31) 
.28|) one 



where A = Adidi. Performing analogous computations starting from the equations (3.26) and (3 
can obtain similar relations for and Cn = qn—iPnfn+i- 

Using the real variables x, y and t 

x = Rezi, y^Im^i, t~Iiez2 (3.32) 
and introducing the real quantities An given by A„ = Re C„ , 

Ai = ^Pn {rn-iqn+1 + Q'ri-lT'n+l) , (3.33) 

one can show that g„ and An satisfy the following system 



idtqn + 2^Qn + ^Anqn = 



kU \qn 



(3.34) 
(3.35) 



where □ — dxx — dyy. This system is nothing other than the DS system jl^ or, to be more precise, the DS-H 
equation, according to the generally accepted classification. 

In such a way we have obtained the following result: solutions of the equations ( p.25| ) - (3.28), belonging 
to the ALH, can be used to obtain solut ions for t he DS equation, i.e. the DS equation can be 'embedded' into 
the ALH. Moreover, when solving ( 3.25| ) - ( 3.28| ), we obtain simultaneously an infinite (for an infinite ALH 
chain) number of solutions for the DS: the pairs ((?„, A„) for all n's solve the DS system. Thus, equations 
( 3.25| ) - ( 3.28| ) can be viewed as sequences of the Backlund transformations. And again, as in the case of 
the 2DTL, we can obtain a wide number of solutions for this (2+l)-dimensional system by, so to say, little 
efforts: without invoking the MIST. But here we come to an important for this approach question. Each 
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solution of the ALH subsystem was shown to satisfy the DS. The question is whether the reverse statement 
is true. The answer is negative. We have seen that all q„'s satisfy the equation of the 0(3,1) cr-model. So, 
solutions that we obtain by 'embedding' into the ALH are only some subclass of all possible solutions of the 
DS. But it turns out that this subclass is rather rich: it contains solitons, 'Wronskian' solutions, finite-gap 
quasiperiodic ones (these solutions are discussed in 0) and many other. It seems to be interesting that 
the pioneering papers by Ablowitz and Ladik which have been written twenty years ago, possess 

everything necessary to construct, say, the N-soliton solutions for the DS equation. 

3.4 The Ablowitz-Ladik hierarchy and the KP equation. 

The result I want present in this section is, in some sense, the most important example of 'embedding' into 
the ALH approach. To derive it one has to consider the s ystem of thr ee first equations of the 'positive' 
subhierarchy, that is the system consisting of the equations ( 2.15| ), ( 2.17 ) and the equation determining the 
third flow, d/dz^: 

diQn = -ip„g„+i (3.36) 

d2qn = -iPnPn+iqn+2 + irn-lPnqnq-n+l + iPnTnqn+i (3.37) 

dsqn = -iPnPn+lPn+2qn+3.+iPnPn+irn+iqn+2 + '^iPnrnPn+iqn+iqn+2 (3.38) 

+irn-lPnqnPn+iqn+2 ~ iPnT^qn+l - ir^-Wnql+l 

+ 2irn-lPnqn+l " ir'^i-lPyiql^qn+l + «f'n-2Pn- iPn^n^n+l 

By very lengthy but straightforward calculations one can verify that for all ri's the quantity 

U = -Kqn-lPnqn+l (3.39) 

satisfies the KP equation: 

di {Adsu + diiiu + I2u dm) = 3922^ (3.40) 
Thus, the KP equation can be 'embedded' into the ALH. 

3.5 The Ablowitz-Ladik hierarchy and the AKNS hierarchy. 

It follows from equations ( 2.15 ) - ( 2.18| ) that the quantities q„ and r„_i satisfy the closed system: 



id2qn + dfqn + 2iq„rn-idiqn = (3.41) 
id2rn-i - dlr^-i + 2ig„r„_iair„_i = (3.42) 

This is one of the forms of the derivative NLSE. To rewrite it in the standard way one can use the gauge 
transform which leads to the following result: the quantities 

Q = g,i exp(-i(?;)), i? = r„_i exp(i(/)) (3.43) 

where 



(izg„(z,Z2)r„_i(z,Z2) (3.44) 



satisfy the system 



idiQ + dlQ + 2idiQ^R = (3.45) 
id2R - dfR + 2idiQR^ (3.46) 

which is the complexified derivative NLSE (I use the term 'complexified' to indicate that Q and R are not 
related by the involution Q = ±^ which is typical for physical applications). Thus, the derivative NLSE can 
be 'embedded' into the ALH. 
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Having obtained this result, one can expect that similar relations exist between the ALH and AKNS 
hierarchy because, using the paraphrase of Kipling's words given by Newell, "The AKNS thing and the 
DNLS string are sisters under the skin" However this cannot be done using only local in g„'s and r„'s 
functions, but is possible in terms of the r-functions of the ALH, 



(3.47) 



(see ( |5.4|) below). It follows from ( 2.15|) - ( 2.18 ) that the quantities Q and R given by 



Tn+l 



Tn-1 



-Tn-l 



(3. 



solve the complexified NLSE 



id2Q + diiQ + 2Q^R = 
-id2R + diiR + 2QR^ = 



(3.49) 
(3.50) 



Moreover, considering also the third flow one can derive that Q and R solve the first of the higher DNLSE's 
as well: 



53Q + 9iiig + 6Qi?aiQ = (3.51) 
d3R + diiiR + 6QRdiR = (3.52) 

This is a demonstration of the fact that the AKNS hierarchy can be 'embedded' into the ALH. 



4 'Embedding' into the 'vectorized' ALH. 

The relations between the 0(3,1) cr-model, 2DTL, DS, KP equations and the ALH, discussed in the previous 
sections were formulated in terms of g„'s and r„'s: solutions of the equations of the ALH yield solutions to 
other problems. Now I want to present some results on the 'embedding' not into the ALH itself but into a 
set of related equations, which are equations for some combinations of the solutions of the linear problems 
(2.4), ( ^ ) associated with the ALH. Such relation is called in literature 'gauge equivalence': NLSE is gauge 
equivalent to the Landau-Lifshitz equation, sine-Gordon — to the Thirring model, DS equation — to the 
Ishimori spin model. 

Consider the matrices cr°, a — 1, 2, 3 defined by 



where (7°(a = 1, 2, 3) are the Pauli matrices 



(4.1) 



1 

1 



-i 

1 



and ^',1 is a 2 X 2 invertible matrix solution of the system {2A), (|2.5| ) 



1 

-1 



(4.2) 



It follows from (|4J|) and (^) that 



(4.3) 
(4.4) 



(4.5) 



The right-hand sides of these equations can be written in terms of the matrices cr" for different a's and n's. 
Moreover, we can obtain the closed system of equations for the matrices (T° with the index a being fixed 
(in what follows we will deal with the matrices cr^) by means of the 'scattering' problem (p3) which can be 



used to relate a'^'s and (Jn±i 
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(4.6) 

(4.7) 



Pn-io-fi-i = (1 + 9„_ir„_i) crfj - 2Ar„_i(T+ + 2A ^g„-if7„ 
where the matrices ct^ are given by 

at^\K±ial) (4.8) 

These identities also enable to express some combinations of the "ALH quantities" Qn, r„ and p„ in terms 
of the "matrix" ones. For example, it follows from (4^) that 



Pn 



1 



1 



2 tr cr„+i(T„ 

For our further purposes it is convenient to use the matrix-vector correspondence 



(4.9) 



5 = ^5aCr° > §={81,32,83) 



(4.10) 



a=l 



where cr" are the Pauli matrices (4.2) 



The equations ( |4.5| ) together with (4.6), (4.7) can be viewed then as a system of DDE's for the matrices 
cr^ or the vectors (?„ which correspond to cr^'s by ( [4.10D , 

erf, <?„ (4.11) 

Namely these equations will play the role of the starting system of DDE's, and from them I will derive 
some interesting difference and differential consequences. Our first example is the 'gauge' analog of the 
0(3,1) (T-model. 



4.1 The Ablowitz-Ladik hierarchy and the Landau-Lifshitz equa- 
tion. 

Using the expressions for V^^, 







A A 2 - r„_i(7„ 



y-l — i ( ^ ^ q-n-lTn ^^n 
Mn-1 



(4.12) 



one can express the derivatives of the matrices cr^ in terms of the matrices given by (4.8) as follows 



(iA/2) dal = 7-„_iCT+ - 9„cr„ , 



(z/2A) Bal 



(4.13) 

(4.14) 



(the symbols d and d stand, remind, for di and di). From these relations, using analogous expressions for 
the derivatives da^, Ba^ and formulae (2.15), ( p. 20 ), one can obtain, after straightforward calculations, 
omitted here, that for every n the matrix 8 = a"^ solves the equation 



308 + I (tr 88 d8) 8 + ^ \X^d8 + X-^d8 , 8] = 
2 2i 



In the vector form equation (4.15) becomes 

8zz + I 8z8z I 8 



X^8z + X'^8z] X 8 



= 



(4.15) 



(4.16) 



which turns out to be the (0-|-2)-dimensional version of the Landau-Lifshitz equation for the isotropic two- 
dimensional classical Heisenberg ferromagnets, 



dt8 = g 8xA8 , 



8"^ = 1 



(4.17) 
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where A is the two-dimensional Laplacian. Indeed, if we restrict ourselves with the stationary structures of 
the form S = S (x — v^t, y — Vyt), then equation (4.17) becomes 
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S xAS 



0. 



Introducing the variables 



z ^ — [x -Vj;t + i{y - Vyt)] , z ^ — [x - vj - i{y ~ Vyt)] 



where v = \v\, the latter can be rewritten as 



S X Szz 



X^S, + X-^Sg = 



(4.18) 



(4.19) 



(4.20) 



with A — exp(i7/2), where the angle 7 is defined by — vcosj, Vy — i;sin7 which is equivalent to (4.16). 
Thus, the (0+2)-dimensional Landau- Lifshitz equation can be 'embedded' into the ALH. The equation 



(4.16) is known to be integrable (its zero-curvature representation one can find in the paper [[18[), and one 
can tackle it by elaborating the corresponding inverse scattering transform, but, as in the cases discussed 
above, we can now obtain a wide range of physically interesting solutions using the already known solutions 
for the ALH. This has been done in the paper [p9|. 



4.2 The Ablowitz-Ladik hierarchy and a multi-cr-field model. 

In this section we will fix the quantity A, 

A = 1 (4.21) 

noting that the case of the arbitrary A's (|A| — 1) can be restored by the rotation of the coordinates x and 
y{z = zi=x + iy) 

Differentiating ( 4.13| ) with respect to z (or (4.14) with respect to z, which leads to the same result) one 
can obtain 



ddan = 2 (2 - p„_i - p„) cr„ + 2 (r„ - r„_i) (t+ -I- 2 (g„_i - g„) (t„ 

Using the identities 

2 [rn - r-n-l) 0-+ + 2 {qn-l - Qn) CTn = PnCr„+i + {pn~l + Pn - 4) (T„ + p„_i(7„_i 



which follow from ( 4^ ) , ( 4/7 ) and ( O ) one can rewrite ( 4.22 ) as 



ddan = P„-lCr„_i - {pn-l + Pn) CTn + Pn(Jn+l 



or, using again (4.13) - (4.14), as 



'^(^<^i + 2 i^^ ^'^fi'^'^D '^n ^ Pn-l<yn-l + {Pn-1 + Pn - 4) O^n + p„(T„+i 



(4.22) 
(4.23) 
(4.24) 

(4.25) 



The last equation turns out to be the field equation of some model which to author's knowledge is new and 
which can be viewed as a multi-field generalization of the well-known 0(3) ct- model. 
The 0(3) a model is described by the Lagrangian 



C = T + C 

where T is the 'kinetic-energy' term in two dimensions 

1 



r=-(^V5,V5j, v = (a,,5,) 

and C takes into account the restriction 

= 1 

Our generalization is, first, in taking instead of one field S an infinite number of fields 



(4.26) 



(4.27) 



(4.28) 
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n = 0,±1, 



i.e., in replacing T (4.27) with 



r 



(4.29) 



(4.30) 



It is easily understood that in absence of interactions between different fields (spins) our generalization would 
be trivial and of little interest both from physical and mathematical standpoints. So, we consider also the 
interaction, the nearest-neighbour one. 



U = 



n+1 



(4.31) 



In what follows our attention will be restricted to the particular pair potential u which has no, so to say, 
physical origin, but is remarkable by the fact that it preserves integrability of the corresponding equations. 
Namely, we will study the potential of the classical Heisenberg ferromagnets model (which is known to be 
integrable, see, e.g., 120, M and to be in close relations with the ALH [p2[): 



-J 



Sn , Sn-^-l 



Summarizing, the model considered is given by 

1 



C = 



n+1 



c 



where 



£' -^A„ (S, 



and A„'s are the Lagrange multipliers, which should be chosen to satisfy the conditions — 1. 
The field equations corresponding to ( 4.33 ) are of the form 



A5„ 



*^ n ~t- ^nSn 



where 



with 



^ n 



dSn 



-U 



2^ ("^"- 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



1 + \ Sn, Sn+l 

and after calculating A„'s they become 

^Sn + {^Sm VSVi^ Sn = -^J |./n-l'5Vi-l + (/n-1 + /n ^ 4) Sn + fnSn+1 1 



(4.37) 



(4.38) 



These equations after identifying /„ and p„ (compare (4.37) and ( [4.9| )) and rescaling the coordinates become 
the vector form of the matrix equations (4.25). Thus, the multi-a -field model described above, which is a 
(2+l)-dimensional system gauge equivalent to the 2DTL, can be 'embedded' into the ALH. 
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4.3 The Ablowitz-Ladik hierarchy and the Ishimori equation. 

The following example is a manifestation of the already known fact that the Ishimori model is gauge equiv- 
alent to the DS equation p^ . 

Using the expressions forl/i±^\ Vt'''^ one can express the derivatives of the matrices cr^ as follows: 



(*A/2) d,al 
(V2A) d,al 

izX/2) d^al 
(V2A) d2al 



(4.39) 
(4.40) 

(4.41) 
(4.42) 



From these relations, using analogous expressions for the derivatives 9jcr^, one can obtain that for every 
n the matrix S = (the index n will be fixed and omitted hereafter) with |A| = 1 satisfies the following 
equations: 



-2id2S +[S,d'^S] = AiwdiS, 



(4.43) 
(4.44) 



where 



W = rn-lQn - A ^, W = Qn-irn - A 



From the other hand, it follows from ( 2.1£ )-( 2.1€ ) that 



idiw = ~idiw = Pn- Pn-l- 



The r.h.s. of ( 4.46 ) can be expressed in terms of the matrix 5 

Pn - Pn^i ^ -Itr {S [diS , BiS]} 
o 



Finally, using the real variables x,y,t (3.32) and the quantity / defined by 



(4.45) 
(4.46) 
(4.47) 



w = -idf (4.48) 
which satisfy the equation A/ = 4 (p„ — Pn-i), equations (4.43), (4.44) and ( 4.46| ) can be presented as 



iSt = -^[S,aS] + iifyS, + USy), 

Af = l:tr {S[S,,Sy]}. 



(4.49) 
(4.50) 



This system, when rewritten in terms of the vector S related to the matrix S by (4.10) becomes the Ishimori 
equation p4] 



St 



Af = (S 



S xas 



Sx X S'l 



(4.51) 
(4.52) 



This is the main result of this section: the Ishimori equation can be 'embedded' into the ALH. 
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5 Functional representation of the Ablowitz-Ladik hi- 
erarchy. 

In the above sections one can find many examples of the relations between the ALH and other integrable 
models. These examples show that the ALH possesses some kind of 'universality'. A hypothesis arises that 
all known integrable models can be 'embedded' into the ALH. So, the main question one should answer 
now is the following: is the ALH in some sense a distinguished hierarchy? It seems to me that there is 
some sense in these words. I cannot give any rigorous proof of such a statement (and even formulate it in 
strict terms). Nevertheless, I think this question is worth studying. The results discussed above are in some 
sense 'empirical' facts: one can easily verify them by simple calculations, but one can hardly find there an 
answer to the question why do such apparently different models turn out to be interrelated. To do this one 
probably has to consider the problem, and hence the ALH, in some more general terms. Namely this was 
the motivation of the work ||] . 

The main result of the paper Q is the representation of the ALH, which has been originally introduced 
as an infinite system of differential-difference equations, as a finite system of difference-functional equations. 
I will not repeat here the calculations outlined there and only write down the main result. The 'positive' 
Ablowitz-Ladik subhicrarchy is shown to be equivalent to following equations: 



[S] 



Tn(z+'- [5]) Pn(z~'- [5]) - Pn (z + [(5]) (z - [6]) 



^[<5]) -a„(z+i[<5])p„ {z-\[5]) 

Tn+l (2 + ^ [5]) Tn-l (z - ^ [(5]) 



(5.1) 
(5.2) 
(5.3) 



Here z stands for (zi, 22, z^, ...), \S\ is the traditional now designation for (5, (5^/2, 6^ /k, ...) , and t„, cr„ 
and pn are r-functions defined by 



In = — , 



I Pn — ^ 



(5.4) 



Expanding (5.1)-(5.3) in the power series in 5 one can obtain all equations of the 'positive' subhierarchy. 
As is seen, in the above formulae only the 'positive' coordinates Zj are used. An analogous representation, 
involving z^'s, can be obtained for the 'negative' subhierarchy. 
Using Hirota's bilinear approach one can rewrite (^.1|)-([5T^) as 



exp 



CTn ■ Tn — Tn • Cr,i — (5 CTn+l • T„_i 
Pn ■ Tn — Tn ■ Pn + S Tn+l ' Pn-l 
Tn ■ Pn — Tn ' T„ + Tn+l ' Tri-l 



^ 



where 



k=l 



and the operators Dzj are defined by 



(- 


d y 


i- 


-A) 




dx' J 


\9y 


dy'J 



(5.5) 



(5.6) 



(5.7) 



- ) u{^,y, ■■■)v{x ,y , ...) 

X —X, y — y,... 

Thus, we have a system of difference-functional equations for an infinite number of triplets of r-functions. 
It has been shown in ^ that from this system one can derive a closed system for only one triplet (ct, p, r) : 
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where cr, p and r stand for (7„, p„ and t„ with rt being fixed and the operator G{5) is defined by 



G{5) = 2 sin 



S Di exp 



(5.8) 



(5.9) 



These equations, when expanded in power series in S yield a hierarchy of partial differential equations of the 
NLSE type. 

This procedure of 'shortening' can be continued by excluding cr and p and obtaining the 'scalar' equations 
for only one function, r, 



H{S) T ■ T = 

where the operator H(S) is given by 

H{d) ^ [2Di - (5 (L»2 + iDii)] exp 



(5.10) 



(5.11) 



Expanding the operator H{d) in powers of S one can obtain an infinite number of operators. The first of 
them is 



4^31 - 3^22 + Dnii 



(5.12) 



i.e. the operator which appears in the bilinear representation of the KP equation. 

Thus, we have a chain: the infinite system (pT|)-( [5^ ) — > the system of three equations ( ^.8[ ) one 
equation (5.10). This chain leads from the ALH, through the NLSE-like hierarchies, to the KP-like one and 
may be useful to understand the place of the ALH among other integrable hierarchies. 



6 Conclusion. 

The results presented in this paper are aimed to demonstrate that the ALH indeed possesses some kind of 
'universality'. One more indication of, in some sense, distinguished character of the ALH, which has not been 
mentioned above and which will be discussed in a separate paper, comes from the theory of the 0-functions. 
It is well known that the ^-functions of the finite genus Riemann surfaces satisfy some functional relations, 
I mean the Fay's trisecant identity There is also a well-elaborated procedure how to derive from the 
latter a number of differential relations satisfied by the 0- functions, and one can easily find in literature a 
lot of identities involving some combinations of some differential operators dj (usually ones that appear in 
the KdV, KP, 2DTL and other equations). It is surprising, but to author's knowledge the following, rather 
natural, question has not been given due attention. Consider the problem of how to express the action of the 
operators dj taken separately on the 0- functions (or some ratios of the functions, 9{C + a) /9{(^ + /3)) without 
invoking other operators dk, {k ^ j), but only in terms of ^-functions, maybe of some other arguments. This 
problem is not very difficult: it can be done using the already known (and widely used) algorithm, by 
expanding the Fay's identity. It turns out that these relations can be written in the form of the equations of 
the ALH. This fact becomes more transparent if we consider this question in the framework of the functional 



representation of the ALH: equations (5.1) - (5.3) in the quasiperiodical case become equivalent to the Fay's 
identity. 

However the word 'universality' often used in this paper should not be understood literally (that is why I 
used it in the quotation marks) because of the following. The relations between the ALH an d ot her equations 



are not, in general, ones of equivalence. The correspondence between the ALH subsystem (|3.l| ) and the field 
equations of the 0(3,1) cr-model (3^) is one-to-one: one can find in j|] the inverse transform (from the field 



equations to the DNLSE-DMKdV chain). At the same time in most of the other examples the situation is 
different. Not all solutions of, say, (2-f l)-dimensional models can be obtained by the 'embedding' into the 
ALH method. And a logical continuation of the works ||, 0, ^ is to study the question of what kind of 
reduction do we have when 'split' some system (2DTL, DS, KP) into few equations from the ALH. In the 
case of, e.g., the DS equation this is easy to do: we can start from the fact that the quantities qn satisfy not 
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only the (2+l)-dimensional DS system (3.34), (3.35) but also the 2-dimensional 0(3,1) cr-model equations, 
and to formulate the ansatz we have implicitly used. However in other cases, the 2DTL for example, this 
problem is somewhat more difficult; I cannot at present write down some, so to say, more simple equations 
which are satisfied by the quantities p„, i.e. 1 cannot formulate, using the language of equations, what 
restrictions do we impose on solutions of the 2DTL when substituting the latter by the ALH equations (3.4). 
I think some progress can be made in the framework of a more general approach, e.g., as one discussed in 
Maybe it will be easier to understand the ALH- (other equations) reductions studying the chain from the 
ALH to the KP equation mentioned at the end of the previous section. To conclude, I want to stress once 
more that this question is not answered yet, while it seems to be one of the most important problems to be 
solved before (if ever) to call the ALH universal (without the quotation marks). 
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